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ABSTRACT 


The  conventional  Optical  Potential  Model  has  been  used 
to  analyze  the  P+4He  elastic  scattering  data  in  the  medium 
energy  range  (  100  MeV  <  Ep  <  1100  MeV  ).  The  model  is 
successful  in  describing  the  data  in  this  energy  range.  We 
find  the  complex  spin-orbit  potential  plays  an  important 
role  in  obtaining  a  good  fit  to  the  data.  An  exchange 
potential  is  also  found  necessary  for  the  reproduction  of 
the  large  angle  data. 

The  behaviour  of  the  volume  integral  of  the  real 
central  potential  indicates  a  logarithmic  energy  dependence. 
A  linear  downward  trend  is  observed  as  the  proton  energy 
increases  and  the  Jr/A  =  0  point  occurs  at  about  450  MeV. 
The  volume  integral  of  the  imaginary  potential  shows  a 
general  upward  behaviour  with  increasing  incident  proton 
energy.  The  r.m.s.  radius  obtained  agrees  well  with  results 
obtained  from  other  sources. 
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CHAPTER  1 


INTRODUCTION 

It  has  long  been  understood  that  to  study  the  Interior 
of  an  atomic  nucleus,  the  wave  length  of  the  projectile  must 
be  comparable  to  or  shorter  than  the  inter-nucleon  spacing 
of  the  target  nucleus.  Therefore,  probing  the  nucleus  with 
high  energy  protons  should  be  more  fruitful  for  nuclear 
structure  research  than  analogous  experiments  at  lower 
energies. 

Recently  a  large  amount  of  effort  has  been  devoted  to 
the  study  of  light  nuclei  with  high  energy  protons,.  The 
pioneering  work  in  this  regard  was  that  by  Palevsky  et  al.  { 
PA67  ) ,  who  obtained  some  interesting  results  for  P+4He 
elastic  scattering  at  1  GeV.  The  differential  cross-section 
was  found  to  exhibit  a  sharp  diffraction  minimum.  Since 
then,  much  experimental  effort  has  been  devoted  to  studing 
the  P+*He  system  in  the  intermediate  energy  range  {  100  MeV 
<  Ep  <  1000  MeV  )  (  PA67,  GO70,  C075,  G068,  VZ75,  BE76, 
MC64 ,  BO 72 ,  SC74,  KL77  )  and  a  large  amount  of  data  were 
reported.  The  most  interesting  features  observed  about  these 
data  are  the  diffraction  minimum  at  about  t  =  -.25(GeV/c)2  ( 
t  is  the  square  of  the  relativistic  four  -  momentum  transfer 
)  for  Ep  >  400  MeV  and  the  rising  back-angle  cross-section 
fcr  proton  energies  below  300  MeV  and  above  600  MeV. 

The  theoretical  analysis  of  the  data  at  medium  energies 
was  performed  exclusively  using  either  the  Glauber 
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Approximation  (  GL59  )  or  Optical  Potential  Model  (  WA57, 

KE59  )  . 

The  Glauber  Approximation,  as  developed  by  Glauber, 
assumes  that  : 

1. )  The  target  nucleons  are  stationary  whilst  the  projectile 
passes  through  the  target  nucleus. 

2. )  Each  collision  between  the  incident  projectile  and  a 
target  nucleon  is  identical  to  the  collision  between  free 
nucleons. 

3. )  The  total  change  in  phase  of  the  particle  wave  function 
is  equal  to  the  sum  of  the  phase  changes  which  occur  at  each 
collision  with  a  target  nucleon. 

4. )  The  Eikonal  Approximation  is  also  assumed  to  be  valid. 


The  scattering  amplitude  F(0)  for  a  spin-0  target  is  given 
by  : 


1 . 1 


where  \  is  the  wave  vector  in  the  center  of  mass  system,  q 

-> 

is  the  momentum  transfer,  b  is  the  impact  parameter  and  x(b) 
is  the  sum  of  the  phases  resulting  from  scattering  with 
individual  nucleons 


j  J 


1.2 


where  s  represents  the 


coordinate 


which 


is 


component  of 
perpendicular 


the  target  nucleon 
to  the  propagation 
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direction  k. 

The  Glauber  Approximation  employed  for  the  early 
analyses  by  Franco  (  FR65  ),  Bessel  and  Wilkin  (  BE67,  68)  , 
and  Czyz  and  Lesniak  ^  CZ67,  68,  70  )  was  able  to  reproduce 
the  dominant  features  of  the  data,  but  failed  to  explain  the 
diffraction  minimum.  Frosch  et  al.  (  FR67  )  pointed  out  that 
the  nuclear  form  factor  should  have  a  zero  at  about  10  fm~2 
due  to  the  result  of  their  analysis.  Hence,  the  gaussian 
form  factor  employed  in  the  early  analyses  is  inadequate  for 
the  *He  nucleus. 

The  Optical  Model  Theory  at  high  energies  was  first 
derived  rigorously  by  Watson  and  collaborators  i  WA51,  53, 
57  )  in  a  series  of  papers  on  multiple  scattering  at  high 
energy.  The  basic  assumption  of  the  theory  requires  the 
Impulse  Approximation  to  be  valid,  which  allows  a 
complicated  many-body  problem  to  be  treated  as  many  two-body 
problems.  He  related  the  multiple  scattering  to  the  Optical 
Model  considering  only  on-shell  elastic  scattering  for  a 
large  number  of  nucleons.  A  briefer  derivation  of  Watson’s 
Model  is  given  in  Chapter  2.  Kerman  et  al.  (K.M.T.)  (  KE59  ) 
took  a  slightly  different  approach,  starting  from  the 
transition  operator  T,  operating  on  the  initial  wave 
function  |  0,  K  >,  where  |  0  >  is  the  target  ground  state 
wave  function.  For  a  properly  antisymmetrized  state  of  the 
target  nucleus 


T  =  Av(  1  +  6T ) 


1.3 
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where  A  is  the  mass  number,  v  is  the  two-body  interaction 
potential  and  G  is  the  Green  Function.  They  defined  an 
effective  two-body  interaction  operator  t  as: 


T  =  V  +  vGt 

to  obtain  the  matrix  element  for  a 
projectile: 


1 . 4 

single  particle 


(A  -  1)<  0 , K 1 1  t 1 0 , i<"  >  1.5 

—y  ~y~ 

for  states  K*  and  K"  of  the  incident  particle. 

The  Optical  Potential  Model  formulated  by  K.M.T.  was 
used  by  Kujawski  \  KU70  )  and  Feshbach  et  al.  (  Fe70,  71,  73 
)  They  included  both  the  first  and  the  second  terms  of  the 
series  in  their  analyses,  where  the  second  term  of  the 
series  should  be  more  sensitive  to  the  nuclear  structure. 
They  were  able  to  obtain  reasonably  good  fits  by  including 
both  spin  and  iso-spin  effects. 

Recently  the  Glauber  Model  regained  some  popularity  due 
to  the  work  by  Wallace  (  WA71,  77  ).  He  expanded  the 
amplitude  in  powers  of  the  difference  between  the  exact  and 
eikonal  propagators.  He  obtained  a  series  having  the  Glauber 
Eikonal  Approximation  as  the  first  term  and  which  converges 
rapidly  to  the  exact  amplitude.  The  so  called  'Wallace 
Corrections'  were  included  in  the  analysis  of  the  1  GeV 
polarization  and  cross-section  data  by  Wallace  \  WA77  )  and 
by  Auger  et  al.  (  AU76  )  to  analyze  the  data  at  0.6,  0.72, 

1.15  and  24  GeV.  These  authors  obtained  excellent  fits  to 
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the  data  by  including  both  spin  and  iso-spin  effects  in 
their  analyses. 

The  Optical  Model  and  the  Glauber  Approximation  can  be 
related  to  each  other  by  means  of  the  Bessel  transform.  The 
idea  was  pursued  by  Dymarz  et  al.  (  DY77  )  who  performed  a 
phase  shift  calculation  where  the  Glauber  profile  function 
was  used  to  generate  an  optical  potential.  These  authors 
were  able  to  obtain  good  fits  to  the  data.  They  compared 
their  approach  with  the  conventional  Glauber  Model  approach 
and  found  the  latter  produced  too  deep  a  minimum. 

Pule  et  al.  {  P075ar  75b,  75c,  76  )  made  an  Effective 
Channel  analysis  of  the  P+4He  elastic  scattering  data.  The 
elastic  scattering  of  high  energy  particles  by  composite 
target  nuclei  is  formulated  in  terms  of  coupled  equations. 
The  effects  of  inelastic  processes  are  represented  by  an 
average  inelastic  channel.  The  theory  incorporates 
approximately  the  effects  of  non-locality,  energy 
dependence,  rescattering  and  absorption  due  to  all  the 
inelastic  channels.  Although  they  ignored  the  spin-dependent 
effects,  fair  agreement  was  still  obtained  for  the 
cross-section  data  at  0.59,  0.72,  and  1.0  GeV 

Phenomenological  Optical  Model  analyses  were  performed 
by  Clark  et  al.  (  AR75,  77  ) .  Using  the  Klein-Gordon 
Equation  or  the  Dirac  Equation,  they  employed  a  simple 
optical  potential  containing  only  a  real  and  imaginary 
central  potential  as  the  fourth  component  of  a  potential 
four  vector.  Even  though  the  spin-orbit  potential  was 


. 
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ignored,  good  fits  to  the  data  were  obtained. 

In  the  lower  energy  region  (  Ep  <  400  MeV  )  most 
analyses  were  done  by  using  either  the  conventional  Optical 
Model  with  the  Schrodinger  Equation  or  the  One  channel 
Resonating  Group  Theory.  The  method  of  resonating  group 
structure  has  only  been  used  to  calculate  cross-sections  for 
scattering  systems  involving  light  nuclei.  It  essentially 
consists  of  building  a  completely  antisymmetrized  wave 
function  for  the  whole  system  (  target  and  projectile  )  from 
channel  functions  which  correspond  to  the  various  ways  of 
grouping  the  nucleons.  This  method  is  a  microscopic  approach 
which  investigates  the  close  interaction  within  the  nucleus 
between  these  groups  as  opposed  to  looking  at  the  nucleus  as 
an  average  potential.  For  Ep  around  100  MeV,  it  is  necessary 
to  introduce  a  phenomenological  imaginary  potential  to 
account  for  all  inelastic  effects,  therefore  the  beauty  of  a 
purely  microscopic  approach  is  lost.  In  general,  the  fits 
obtained  by  the  Resonating  Group  Theory  {  RA74  )  and  by  the 
Optical  Potential  Model  are  equal  in  quality.  However,  the 
Optical  Potential  Model  is  more  convenient  for  analyses  in 
the  same  energy  range. 

The  question  of  the  significance  of  the  spin-orbit 
potential  has  been  raised  in  many  of  the  analyses  mentioned 
above.  Clark  et  al.  (  AR76  )  have  analyzed  the  data  at  a 
number  of  energies  above  500  MeV.  There  has,  as  yet,  been  no 
attempt  to  perform  an  energy  dependent  analysis  of  all  the 
experimental  data  from  100  to  1000  MeV.  In  this  thesis  we 
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want  to  examine  the  validity  of  using  a  conventional  Optical 
Potential  Model  to  fit  the  P+4He  elastic  scattering  data  and 
the  effects  of  the  spin-orbit  potential  in  the  medium  energy 
range. 

In  chapter  2,  we  give  a  brief  derivation  of  Watson*s 
multiple  scattering  model  and  some  of  its  characteristics. 
In  chapter  3 ,  we  present,  a  complete  description  of  the 
method  of  analysis.  In  chapter  4,  we  compare  the 
experimental  data  with  our  results  at  each  corresponding 
energy.  In  chapter  5f  we  discuss  our  results  as  compared  to 
other  methods  of  analysis.  The  computing  program  •SEEK*  is 
described  in  Appendix  1. 


•- 
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CHAPTER  2 


THEORY 

2, 1  Introduction 

The  Optical  Potential  Model  has  been  used  extensively 
to  analyze  elastic  scattering  experiments  and  direct 
reaction  experiments.  The  theory  is  successful  in  describing 
both  high  and  low  energy  nuclear  experiments.  The 
construction  of  the  model  is  purely  phenomenological  and 
based  on  the  assumption  that  the  atomic  nucleus  can  be 
treated  as  one  single  potential  instead  of  a  many-body 
interaction  between  the  incident  particle  and  the  target 
nucleons.  The  question  about  the  physical  meaning  of  such  a 
phenomenological  model  has  been  asked.  Such  a  profound 
question  has  been  answered  by  Watson  and  collaborators  who 
provided  the  theoretical  justification  for  the  use  of  the 
Optical  Model.  In  a  series  of  papers  on  the  theory  of 
multiple  scattering,  Watson  and  his  collaborators  indicated 
the  manner  in  which  the  Optical  Potential  for  high  energy 
nucleon-nucleus  scattering  could  be  calculated  presupposing 
a  knowledge  of  nucleon-nucleon  interactions.  Since  we  are 
only  interested  in  elastic  scattering,  we  shall  give  a 
simple  derivation  of  the  Optical  Model  for  high  energy 
elastic  scattering.  The  derivation  follows  closely  the  work 
by  Watson  et  al.  (  WA51,  57  )  and  by  Jones  (  J063  ). 
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2.2  Theory  of  the  Optical  Model  Potential 

Let  H.  be  the  Hamiltonian  of  the  target  nucleus,  and  e 
z  n 

be  the  energy  eigenvalue.  The  equation  of  motion  for  the 

nucleus  is  : 


2.  1 


The  target  Hamiltonian  is  assumed  to  be  hermitian  and  the 
states  |  £  >  are  completely  antisymmetric  and  form  a 
complete  orthogonal  set  with  ground  state  j£Q>.  The  equation 
of  motion  for  the  incident  particle,  when  there  is  no 
interaction  with  the  nucleus,  is 


H  |T< ’>  =  E  lie '> 
P1  P1 


2.2 


where  H  is  the  free  projectile  Hamiltonian,  E  is  the 
P  P 

kinetic  energy  and  is  the  momentum  of  the  incident 

particle. 

The  first  essential  approximation  is  to  neglect  any 
antisymmet rization  between  the  incident  nucleon  and  the 

target  nucleons.  If  thefe  is  no  antisymmetrization  with  the 
incident  particle,  one  can  assume  that  the  origin  of  the 
Optical  Model  Potential  should  coincide  with  the  center  of 
mass  of  the  target  nucleus  for  any  mass  number  A.  This 

assumption  is  common  to  all  phenomenological  Optical  Model 
analyses.  However,  if  the  incident  nucleon  is 

antisymmetrized  with  the  target  nucleons,  it  is  reasonable 
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to  require  that  all  A+1  nucleons  should  have  a  common 
single-particle  potential  with  a  specified  origin.  This 
condition  does  not  seem  to  be  possible  due  to  the  recoil  of 
the  target  nucleons.  Watson  et  al.  (WA57  )  investigated  the 
effects  of  antisy mmetrization  and  found  those  effects  to  be 
small  and  unimportant. 

The  Schrodinger  Equation  we  have  to  solve  is  : 

(H0  +  V)|T>  =  E|T>  2, 3 

where  H  =  H.  +  H 
0  f  p 

E  =  e  +  E 

n  p 

V  is  the  interaction  potential  between  the  projectile  and 
the  target,  and  j  Y >  =  |  £s  £'  >.  Because  both  |£  >  and 

I  >  are  complete  sets  their  product  states  |  ¥  >  must  be 

both  complete  and  orthogonal.  The  equivalent  form  of  | V  >  is 


y>  =  |£  ,Tc’>  +  G(E)V | ¥> 
o 


2.4a 


where 

G(E)  =  (E-H0+1e)_1 


2.  4b 


The  second  term  in  [2. 4a)  includes  the  boundary  condition  of 
out-going  scattered  waves  only,  and  |  £  , k 1  >  is  the  initial 
state  wave  function,  consisting  of  a  plane  wave  and  the 
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ground  state  of  the  target  nucleus. 

The  transition  probability  of  the  system  depends  on  the 
matrix  element  of  T,  which  is  given  by: 


T  =  <i  ,1c '  I V I  Y> 

o 


2.5 


For  elastic  scattering,  we  can  separate  the  ground  state 
component  by  a  projection  operator  3  such  that. 


$» 


j _ r 

(2ir)3J 


die"  I?  ,£"><C  ,1c" 

0  0 


=  3  V> 


2.6 


The  Optical  Model  operator  is  the  interaction  operator  which 
leads  to  the  elastic  scattering  state  vector  |  $  >  only.  The 
Optical  Model  Operator  is  defined  by  the  equation 


<D>  =  |£  ,P>  +  G(E)U|<f» 


2.7 


The  operator  U  must  be  such  that  it  is  incapable  of  changing 
the  target  ground  state  wave  function  into  other  excited 
states.  Therefore  we  can  write 


<C,lc|U|  5  ,£'>  =  6,  ,  <llU|l<'> 

0  s  t, 


2.8 


and  U  is  referred  to  as  a  single  particle  operator. 


■ 


■ 
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Following  the  derivation  by  Watson,  we  define  an 
operator  F  such  that 


y>  =  F  U» 


2.9 


Then 

F  =  1  +  G(E)  (V-U)F  2.10 

where 

G(E)  =  (  E  -  U  -  Hq  +  ie  )_1  .  2.11 

From  (2.5),  the  elastic  scattering  amplitude  can  be  written 
as 


<£  I U I $>  =  <£  |  V I  ¥> 

o  o 

=  <£  Jc'IVFl^ 

o 

Substituting  for  F  from  (2.10),  we  have 

<c  ,f'|V[l+G(V-U)+G(V-U)G(V-U)+  •  •  -]|l» 
0 


2.  12 


2.  13 


2.  14 


From  the  equation  of  G  (E) ,  (2.4b)  can  be  written  as  : 


G(E)  =  [  G”1 (E)  +  U  I'1 


2.  15 


Using  matrix  algebra  one  can  rewrite  (2.15)  as 


' 


* 
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G  =  G  -  GUG  +  GUGUG  -  •  •  • 
or 


2.  16 


1  -  GU  =  1  -  GU  +  GUGU  -  GUGUGU  +  •  •  • 


Now,  rearrange  the  series  (2.14)  as  follows  : 

V  (  G  -  GUG  +  •  •  •  ) V  •  •  • 

•  •  •  (  G  -  GUG  +  •  •  •  )V(  1  -  GU  +  •  •  •). 
Using  (2,17)  and  (2.18)  we  can  rewrite  (2.14)  as 


2.  18 


<U>  =  <  (  V  +  VGV  +  •  •  •  )(  1  -  GU)  >  # 


2.19 


The  intermediate  states  which  may  occur  in  the 
scattering  process  are  such  that  either  the  incident 
particle  or  the  target,  or  both,  may  be  in  excited  states. 
It  is  convenient  to  separate  the  states  in  which  only  the 
incident  nucleon  is  excited,  with  the  target  remaining  in 
the  ground  state.  This  can  be  done  by  a  projection  operator 
defined  by: 


a  =  1  - 


% 


5  ,1c" 

0  0 


2.20a 


and  is  related  to  3  by* 


1  =  a  +  8 


2.20b 


' 


•  • 
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Now,  the  equation  (2.14)  can  be  written  as: 
<U>  =  < | Y [  1  -  G  ( a+  3 )  U  ]|> 

=  < | Y[  1  -  3GU  ]|> 


2.21 


Recall  that  U  only  acts  on  the  target  ground  state  and  a 
direct  iteration  of  (2.21)  gives 


<U>  =  <|  Y  -  YBGY  +  Y6GY3GY  +•••!> 


2.  22 


where 

Y  =  V  +  VGV  +  •  •  • 

=  V  +  V(a+3)GV  +  •  •  • 
Regrouping  Y  gives 

Y  =  L  +  L3GL  +  L3GL3GL  +  •  •  • 

where 


2.23 


2.24 


2.25 


L  =  V  +  VaGL  +  •  •  • 


2.26 


Equation  (2.26)  shows  that  the  intermediate  states  are  all 
such  that  the  target  and  incident  particles  are  excited. 
However,  in  the  series  (2.25)  all  the  intermediate  states  of 
the  target  nucleus  are  ground  states,  only  the  incident 
nucleon  being  excited. 

The  potential  V  is  taken  to  be  the  sum  of  two-body 


. 


■ 
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interaction  potentials,  namely 


A 


2.27 


At  this  stage  we  introduce  a  two-body  operator  t. : 


t .  =  v.  +  v,aGv,  +  •  • 


2.28 


v.  +  v-aGt. 


2.29 


The  operator  t.  represents  the  scattering  of  the  incident 
particle  by  the  i*'*1  target  nucleon.  This  operator  is  not 


equal  to  the  "free"  transition  operator  T  in  the  usual  sense 
due  to  the  presence  of  aG  factor.  However,  at  high  energies 
the  impulse  approximation  is  assumed  to  be  valid,  hence 
these  two  exprssions  become  the  same. 

Now  substitute  (2.25)  into  (2.22)  after  a  direct 
expansion.  We  obtain 


<U>  =  <L> 


2.30a 


or 


<U>  =  < 


> 


2.30b 


In  the  series  (2.30)  we  can  see  the  Optical  Potential  as 
being  represented  by  a  series  of  two-body  interaction 


. 

• 
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potentials •  The  first  term  of  the  series  is 


2.31 


where  k  is  the  initial  and  ic  *  is  the  final  momentum  in  the 
center  of  mass  system  of  particle  and  nucleus.  If  we  assume 
the  validity  of  the  Impulse  Approximation  and  complete 
an t i-symmetrization  of  the  ground  state  wave  function  of  the 
target  nucleus,  then  all  the  terms  in  the  sum  are  equal  and 
(2.31)  can  be  written  as 


<£'  I  It  t,- 1£>  =  A<i< 1 1 1  | t> 

i  1 1  1  o 1 


A  f 


w 


dP  dP'  dj 

0 


o 


2.32 


where  P  and  P*  are  the  initial  and  final  values  of  the 
momentum  co-ordinate  of  the  target  nucleon.  The  symbol 
denotes  the  remaining  A-1  momentum  co-ordinates  which  we 
assume  to  be  unchanged  by  the  collision.  Before  we  proceed 
further,  we  should  point  out  that  the  essential 

approximation  done  here  is  to  replace  tj  ,  the  effective 
two-body  operator,  by  ,  the  transition  operator  for  the 
scattering  of  two-free  nucleons. 


. 


. 
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For  high  energies  the  dependence  on  the  motion  of  the 
target  nucleon  can  be  neglected  and  the  matrix  element  of  tQ 
is  given  by: 


<#',£' |t0|?t>  =  <£■ |t 


2.33 


Substituting  (2.33)  into  (2.32),  we  have 


A 


(2tt) 


?A 


d y  d?  V*  (y,P-q)<P' |  tJP>4^(y,P) 


2.34 


where  §  =  k*-  -  k,  is  the  momentum  transfer.  The  Fourier 

Transform  of  the  target  wave  function  )  to  momentum 

o 

space  is 


V5(y,?) 

0 


->  -i  (P*r  '+y*n 1 ) 

dn  dr  i(n  ,r  )  e 


2.35 


and  implies  that  i2. 32)  can  be  written  as:- 


-> 


<P'|U|P>  =  A<1< '  I tQ 1 1<>  dn'  dr'  e~1C|’r  l^n  1  ,r' )  | 


A <t'  |  tQ|  1<> 


->■ 


dr'  p(r')  e"iq*r 


2.36 


where 


dn'  |^(n'  ,r')  |2 
o 


2.37 


. 


' 
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and  p  is  the  onebody  density  distribution  normalized  by:- 


f 

. 


p(r')  dr 


l 


=  1 


2.38 


If  U  is  local  and  energy  dependent,  then 


<£'|U|£>  =  f  dr'  U(k,r' ) 


2.39 


and  for  g  =  0 


j  dr'  U(k,r')  =  A<1< 1 1  tQ | lt>  #  2.40 

Hence,  we  can  see  the  volume  integral  of  the  optical 
potential  is  simply  related  to  the  nucleon-nucleon 
scattering  amplitude. 

To  understand  the  significance  of  the  imaginary 
potential,  consider  q  =  0;  the  forward  scattering  amplitude 
is  given  by 

f(0)  -^<*|t0|t>  2.41 

4tt 

By  the  Optical  theorem,  we  have 

Imf(O)  =  — |f|  .  2.42 

8tt 

where  V  is  the  relative  velocity  of  the  colliding  nucleons 
and  is  the  effective  total  nucleon- nucleon  cross-section 


' 

..  ■ 

■ 


■ 
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averaged  over  neutrons  and  protons.  Equating  (2.42)  and 
(2.41) 


Im(-^-<l(|tJl<>)  =  -22_|V 

4ir  0  8tt 


2.  43 


or 


Im(  }  = 

Ap(r) 


which  implies 


2.  44 


-  im  u(k.r')  -  JSsMH 

2 


2.45 


A|tf 

2X 


2.46 


where  p  is  the  nucleon  density  for  the  target  nucleus.  We 
can  see  that  the  imaginary  potential  is  simply  related  to 
the  mean  free  path  X  inside  the  nuclear  matter.  However,  one 
should  note  that  the  Exclusion  Principle  has  not  beer- 
accounted  for  here.  A  detailed  treatment  has  been  given  by 
Lane  et  al.  (  LA65  ) . 


■ 


CHAPTER  3 


METHOD  OF  ANALYSIS 
3. 1  Optical  Potential  Model 

The  P  +  *He  elastic  scattering  data  in  the  medium  energy 
range  (  100  MeV  <  Ep  <  1000  MeV  )  are  analyzed  using  the 
Optical  Potential  Model.  The  target  nucleus  is  assumed  to  be 
spherical  and  to  have  spin-0.  This  Optical  Model  Potential 
is  an  approximation  to  the  true  potential  which  represents 
the  Optical  Model  Operator  U  defined  in  Chapter  2.  The 
Optical  Model  Potential  used  here  has  the  form: 


V  .  (r)  =  V  (r)  -  V  f  (r)  -  iW  f  (r) 
opt  C  0  0  V  V 


P  df. 

+  -(  V  +  iW  )  , 
r  s*o  s*o  dr 


s  •  o  rt  -± 


a-o) 


3.  1 


where  V  (r)  denotes  the  coulomb  potential  due  to  a  uniformly 
c 

charged  sphere  of  radius  R  and  charge  Ze: 


Vc(r)  - 


Ze 


2R  (  3_(R  )2  5 
c  c 


r  <  R 
—  c 


= 

r 

where  the  radius  Rr  is  taken  to  be  1.36A1/3  fm  which  is 

V/ 

taken  to  be  constant  throughout  the  analysis.  The  quantity 
VQ  denotes  the  real  central  potential  and  Wv  denotes  the 


3.2 


r  >  R. 
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volume  part  of  the  central  imaginary  potential.  The 
spin-orbit  potential  is  of  the  usual  Thomas  form  with  Vs ^ 
as  the  real  part  and  as  the  imaginary  part.  The  shape 
of  the  potential  is  expressed  through  a  function  f . (r)  of 
the  Woods-Saxon  form: 


f^r)  =  [  1  +  exp(-—  ]  -1  3.3 

ai 

In  general,  the  radius  and  diffusness  of  the  form 
factor  are  different  for  different  parts  of  the  potential. 
During  the  course  of  our  search,  we  found  no  evidence  of  a 


surface 

peaked  imaginary 

term. 

consequently 

it 

is  not 

included 

in  equation  (3.1) 

and  is 

not  used  in 

the 

latter 

part  of  the  analysis. 

3.2  Data. 

There  are  a  large  number  of  P+4He  elastic  scattering 
data  available  for  cross-section.  However,  some  of  these 
experimental  results  have  changed  very  rapidly  during  the 
last  few  years.  In  order  to  obtain  any  meaningful  results, 
we  have  to  discriminate  against  those  data  that  are  known  to 
be  inconsistent  with  other  authors  work.  In  table  1  we  give 
a  list  of  experimental  data  selected  for  this  search. 


. 


*4 
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3.3,  Relativistic  Corrections 

When  using  the  non-relativistic  Schrodinger  equation  to 
analyze  medium  energy  datar  it  is  necessary  to  correct  for 
relativistic  effects.  If  these  effects  were  not  corrected, 
the  analysis  would  give  misleading  information  about  the 
potential  depth  and  the  geometry.  The  problem  of 
relativistic  corrections  could  have  been  avoided  if  the 
Klein-Gordon  equation  or  the  Dirac  equation  were  used 
instead  of  the  Schrodinger  Equation.  Even  so,  problems  still 
exist.  In  the  first  case,  if  one  were  to  use  the 
Klein-Gordon  equation  to  describe  the  scattering  of  spin-1/2 
particles  scattering  off  a  spin-0  target,  then  one  would 
have  to  ignore  the  spin-orbit  effects.  The  effects  of  the 
spin-orbit  potential  are  clearly  essential  for  the 
reproduction  of  polarization  data  as  indicated  in  the  lower 
energy  analyses.  In  particular  many  authors  have  pointed  out 
that  the  spin-orbit  force  is  essential  in  the  analyses  of 
F+* He  scattering  data.  A  more  accurate  way  to  include  both 
spin  and  relativistic  effects  is  by  using  the  Dirac 
equation.  The  Dirac  equation  with  an  explicit  spin-orbit 
potential  has  not  yet  been  solved.  There  are  also 
ambiguities  in  the  definition  of  the  Optical  Potential  that 
have  not  yet  been  fully  explored. 

Elton  i  EL66,  IN74  )  pointed  out  that  one  may  correct 
for  the  relativistic  effects  by  calculating  the  wave  number 
«K*  in  the  center  of  mass  frame  using  relativistic 


. 


. 


' 

r , .  .  •  ■  ■  ■  *  • ' 


■ 


23 


kinematics.  Such  a  wave  number  is  given  by 


K  =  ( 


Tlab(  Tlab+2mi) 

(  m-|+m2  )2  +  2m2T 


1/2 


lab 


3.4 


where  m  and  m  are  the  rest  masses  of  the  projectile  and 
target  respectively  and  is  the  kinetic  energy  of  the 
incident  projectile  in  the  laboratory  frame.  One  then  works 
backward  using  non-relativistic  mechanics  to  obtain  an 
incident  energy  to  use  as  a  parameter  in  the  automatic 
search  code  'SEEK*  (  SE65  ). 

3.4.  Search  Procedure. 

In  the  present  analysis  of  P+4Re  elastic  scattering 
data,  the  final  parameters  were  determined  generally  by 
minimum  values  of  the  total  chi-square.  This  quantity  is 
given  by: 

X2  =  +  X*  3.5 


where  y2 


is  defined  as: 


Na 

I  r- 

i=l 


ex 


(e,)  - 


°th^0i^ 


toex{01) 


3.6 


and  cr  ( ©1 )  ,  °th(0i)  *  ana  A<WV  are  respectively,  the 
experimental,  theoretical  and  experimental  error  values  for 


. 


. 


■ 


' 
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the  cross-section  at  a  center  of  mass  angle  0.  .  A  similar 

i 

expression  holds  for  the  x2  ,  the  x2  associated  with  the 
polarization  data. 

In  the  beginning  of  the  search,  some  criteria  were  set 
to  eliminate  some  of  the  ambiguities  that  might  arise  from 
considering  a  pure  minimum  y2  fit  to  the  experimental  data 
alone.  The  conditions  are  : 

1. )  The  calculated  reaction  cross-section  must  agree 
reasonably  well  with  the  experimental  data  whenever  it  is 
available  (  IG67  ). 

2. )  Parameters  must  exhibit  a  reasonable  continuity  as  a 
function  of  the  incident  proton  energy. 

3. )  Whenever  possible,  cross-section  data  and  polarization 
data  are  analayzed  simultaneously. 

The  first  two  conditions  are  used  to  eliminate  any 
irregular  values  of  the  optical  model  parameters  in  the 
search  processes.  Therefore  some  fits  with  smaller  value  of 


were 

discarded. 

In 

phenomenological  Optical 

Model 

analysis,  an 

ideal 

of 

data  should  contain 

bot  h 

polarization 

and 

cross-section  data  for  simultanteous  analysis.  The  presence 
cf  polarization  data  helps  to  regulate  the  variations  of  the 
spin-orbit  term.  The  desire  of  having  both  polarization  and 
cross-section  data  at  the  same  energy  can  not  be  satisfied 
by  all  data  sets.  Therefore  to  control  the  variations  in 
spin-orbit  potential,  we  included  polarization  data  from  a 
neighbouring  energy  with  existing  cross-section  data  into 


.  a  ♦  -  •:  ac.  t  l'.  *  ■ 

■ 


■ 


25 


the  analysis  and  vice-versa. 

The  optical  model  parameters  at  156  MeV  found  by 

Comparat  et  al. (  C075  )  were  used  as  starting  parameters  of 

the  analysis.  The  parameters  were  improved  by  the  computer 

code  ’SEEK*  at  156  MeV.  Attempts  were  made  to  maintain  the 

improved  geometry  set  at  a  constant  value  for  different 

proton  energies.  This  procedure  was  carried  out  until 

reasonable  fits  to  the  experimental  data  could  not  be 

obtained  by  varying  the  potential  parameters  alone.  In 

keeping  the  geometry  constant  we  reduced  the  number  of 

variables  from  eleven  to  four  (V  ,  W  ,  V  and  W  ) . 

o  v  s  •  o  s  •  o 

We  hope  this  kind  of  practice  will  provide  a  better  picture 
of  energy  dependence  of  Optical  Potential  parameters. 

The  improved  geometry  set  obtained  from  the  above 
procedure  failed  to  reproduce  the  P+*He  data  above  312  MeV. 


A  new  set 

of  geometry 

was  developed 

by 

choosing  an 

appropriate 

set  of  a0 

and  3^  to  match 

the 

r.m.s.  radius 

obtained  from  e+*He  elastic  scattering  data.  Final 
variations  were  made  on  aQ  rr0  and  other  parameters  to 
minimize  x2  •  This  new  set  of  geometry  was  used  for  proton 
energies  from  500  MeV  to  1100  MeV. 

For  large  scattering  angles,  the  Optical  Model 
Potential  (3.1)  is  not  able  to  reproduce  the  experimental 
data.  In  an  attempt  to  remedy  this  deficiency,  a  Majorana 
Exchange  Potential  term  (  GF72  )  is  added  to  the  potential. 


* 


. 


, 
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The  Exchange  Potential  is  defined  as: 


3.  7 


Now  the  modified  optical  potential  model  is  given  as: 


U(r)  -  VQpt(r)  ♦  (-D\xfex(r)  .  3.8 

Initially  Optical  Model  searches  on  large  angle  data 
were  performed  at  156  MeV  by  starting  with  the  best  fit 
parameters  from  the  forward  angle  fit(  Oc.m.  <  100°  ), 
allowing  only  the  parameters  associated  with  the  exchange 
potential  to  vary  in  order  to  fit  the  data  including  the 
back-angles.  After  convergence  of  the  x2#  final  variations 
on  all  potential  parameters  were  allowed  to  minimize  the  x2« 
In  general,  the  final  changes  were  small.  The  exchange 
potential  geometry  obtained  is  generally  kept  constant  for 
further  analyses.  The  geometry  of  the  Exchange  Potential  was 
allowed  to  vary  only  if  the  geometries  of  V  t  (r)  were 


allowed  to  change 


- 


■ 
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Table  1 


Data  selection  for  the  P+*He  Optical 
Model  analysis  in  the  energy  range  100  -  1150 


Tp  Types  of  Reference 

(  MeV)  data 


100 

147 

156 

203 

312 

350 

560 

650 

650 

720 

720 

800 

1030 


cross-section 

l  GO 70  ) 

polarization 

(  C059  ) 

cross-section 

(  CO  7  5  ) 

polarization 

(  G059  ) 

z 

cross-section 

polarization 

(  CH56  ) 

z 

cross-section 

cross-section 

(  SC  76  ) 

cross-section 

{  KL7 7  ) 

z 

polarization 

cross-section 

l  SC  76  ) 

cross-section 

l  BE76  ) 

cross-section 

(  VE75  ) 

cross-section 

l  MC65  ) 

z 

polarization 

cross-section 

(  KL  7  7 

z 

polarization 

cross-section 

(  KL77 

z 

polarization 

cross-section 

(  SC  76 

MeV 


1150 


- 


■ 
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CHAPTER  4 


RESULTS 

The  P+4He  data  have  been  analyzed  by  the  Optical 
Potential  Model.  The  automatic  search  for  minimum  x2  was 
performed  by  the  computer  code  ‘SEEK*  (  App.  1  ),  which 
allows  automatic  search  on  all  parameters.  Since  we  have 
fixed  the  geometry  parameters  in  most  cases,  the  number  of 
search  parameters  are  signif icantly  reduced  and  the  cost  of 
the  automatic  searches  is  also  reduced.  The  quality  of  our 
fits  is  based  on  the  over  all  fit  to  the  data  set.  The  total 
X2  is  only  used  as  guidance  when  performing  automatic 
searches.  The  results  obtained  from  the  searches  taking  into 
account  the  criteria  discussed  in  Chapter  3  are  detailed  in 
tables  2  and  3.  In  figures  1  to  16,  the  theoretical  curves 
calculated  from  the  best  fit  parameters  together  with  the 
experimential  data  are  presented  at  each  individual  energy. 

1 .)  100  MeV  DATA  <  Fig.  1  ) 

The  data  set  is  taken  from  Goldstein  et  al.  (  GO70  )  •  It 

contains  only  the  differential  cross-sections.  The  starting 

search  parameters  were  taken  from  the  156  MeV  best  fit 

result.  The  method  used  in  the  search  was  to  consider 

related  parameters  in  pairs  i  i.e.  V  and  W  and,  V  and 
r  o  v  s*o 

Ws.0  pairs) ,  until  the  Xz  converged.  The  Optical  Model 
calculation  gives  a  good  fit  to  the  forward  angle  data,  but 
an  exchange  potential  is  required  to  reproduce  the  large 


28 
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angle  data.  The  Optical  Model  calculation,  including  the 
exchange  effects,  is  shown  in  fig.  1. 

2. )  156  MeV  DATA (  Figs.  2  and  3  ) 

The  147  MeV  polarization  data  i  C059  )  and  the  156  MeV 
cross-section  data  (  C075  )  were  considered  together  in  the 
analysis.  An  exchange  term  is  included  for  the  reproduction 
cf  large-angle  data  (fig.  2).  The  polarization  data  at  the 
forward  angles  are  well  describe  by  the  calculation,  but  the 
calculation  can  not  reproduce  the  large  angle  data  with  the 
same  quality  (  fig.  3  ) . 

3. )  203  MeV  DATA  (  Fig.  u  ) 

The  data  set  is  obtained  from  Gotow  et  al.  (  G059  )  .  It 
consists  of  both  cross-section  and  polarization  data.  The 
starting  parameters  are  taken  from  the  156  MeV  best  fit 
result.  Using  the  method  described  above  we  find  the 
calculation  gives  a  good  representation  for  both 
cross-section  and  polarization  data  (fig.  4) . 

4. )  312  MeV  DATA  (  Fig.  5  ) 

The  data  set  is  obtained  from  Chamberlain  et  al. 

(  CH 56  )  It  contains  both  cross-section  and  polarization 
data.  The  starting  parameters  are  taken  from  203  MeV  best 
fit  result.  The  calculation  provides  a  good  fit  to  the  data. 
We  should  note  that  for  proton  energies  at  100,  203  and 
312  MeV,  the  number  of  variables  were  limited  to  four  only 
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and  the  geometry  set  is  the  same. 

5. )  350  MeV  DATA  (  Fig.  6  ) 

The  data  set  is  that  of  the  Saclay  Group  (  SC76  )  .  It 
contains  both  cross-section  and  polarization  data.  The 
geometry  used  in  the  preceeding  analyses  failed  to  provide  a 
good  fit  to  the  data.  A  new  set  of  geometry  parameters 
(table  2)  were  introduced  into  the  search.  It  was  then 
possible  to  get  a  good  fit  to  the  data. 

6. )  560  MeV  DATA  (  Figs.  7  and  8  ) 

The  560  MeV  data  set  (  KL77  )  contains  both 
cross-section  and  polarization  data.  Using  the  method 
described  before,  we  find  the  calculation  gives  good  fits  to 
both  the  cross-section  data  and  up  to  the  first  minimum  (  t 
=  - . 2  (GeV /C ) 2  )  for  polarization  data. 

7. )  650  MeV  DATA  (  Fig.  9  ) 

The  data  set  includes  two  separate  sets  of  cross-section 
data.  The  forward  angle  data  is  from  Aslanides  et  al. 
<  SC76  )  and  the  back-angle  data  is  from  Berger  et  al.  ( 
BE76  )  .  The  calculation  gives  fair  agreement  with  the  data 
by  including  the  exchange  effects  (  fig.  9  )  . 

8. )  720  MeV  DATA  (  Figs. 10  and  11  ) 

The  data  set  is  made  up  of  two  cross-section  data  sets 
and  a  set  of  polarization  data  (  VE75,  MC65  ).  Employing  the 
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usual  method,  we  find  the  calculation  gives  an  excellent  fit 
to  the  cross-section  data.  A  fair  agreement  is  also  obtained 
for  the  experimental  polarization  data. 

9. )  800  MeV  DATA  (  Figs.  12  and  13  ) 

The  data  set  includes  both  cross-section  and 
polarization  data  (  KL77  ) .  By  using  the  best  fit  parameters 
from  720  MeV  as  starting  parameters  we  were  able  to  obtain 
a  good  fit  to  the  cross-section  data  and  a  fair  agreement 
for  the  polarization  data. 

10. )  1030  MeV  DATA  (  Figs.  14  and  15  ) 

The  data  set  includes  both  cross-section  and 
polarization  data  (  KL77  )  .  Taking  the  800  MeV  best  fit 
results  together  with  the  usual  search  method,  we  find  the 
calculation  provides  good  fit  to  the  cross-section  and  small 
angle  polarization  data  (  t  <  -. 2(GeV/C) 2  ). 

11. )  1150  MeV  DATA  (  Fig.  16  ) 

The  data  set  contains  cross-section  data  (  SC76  )  only. 
By  the  usual  method  of  calculation  excellent  fit  to  the  data 
is  obtained.  For  proton  energies  at  560,  650,  720,  800,  1030 
and  1150  MeV  the  search  parameters  were  limited  to  four  and 
an  extra  degree  of  freedom  (Exchange  Potential)  was  added 
for  the  fitting  of  large  angle  data. 
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Fig  1.  Optical  Model  fit  to  100  MeV  data. 
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Fig  2.  Optical  Model  fit  to  156  MeV  cross-section  data. 
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.  Optical  Model  fit  to  147  MeV  polarization  data 
The  theoretical  curve  is  calculated  for  Ep=l 56  MeV 
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Fig  4.  Optical  Model  fit  to  203  MeV  data. 
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Fig  5.  Optical  Model  fit  to  312  MeV  data. 
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Fig  6.  Optical  Model  fit  to  350  MeV  data. 
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Fig  7.  Optical  Model  fit  to  560  MeV  cross-section  data. 
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Fig  8.  Optical  Model  fit  to  560  MeV  polarization  data. 
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Fig  9.  Optical  Model  fit  to  650  MeV  data. 
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Fig  10.  Optical  Model  fit  to  720  MeV  cross-section  data. 
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Fig  11.  Optical  Model  fit  to  720  MeV  polarization  data. 
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Fig  12  Optical  Model  fit  to  800  MeV  cross-section  data. 
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Fig  13.  Optical  Model  fit  to  800  MeV  polarization  data. 
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Fig  14.  Optical  Model  fit  to  1030  MeV  crOss-section  data. 
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15.  Optical  Model  fit  to  1030  MeV  polarization  data. 
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Fig  16.  Optical  Model  fit  to  1150  MeV  data. 
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Optical  Model  Geometry 


For 

Fp  = 

100 

to  312  MeV 

a 

0 

=0.  68 

f  m 

a. =0.317 

1 

fm  a  =.265fm 

s-o 

a  =.408fm 
ex 

ro 

=  1.56 

f  m 

r -j  =  1 . 540 

fm  r  =.802fm 

S  .0 

reX=1.5«£m 

rc 

=  1. 36 

f  m 

IT  •  Id  •  S  •  — 

3.073  fm 

Fcr  Fp=350  MeV 


aQ  =0.  322f m 

a-j=0.305fm  as#0  =.  226fm 

rQ  =1 . 064f m 

r^  =  1.111fm  rs.o=0.831fm 

rr  =1 . 3  6  f  m 

r.m.s.  =  1.77  3  fm 

For  Fp=  440  to 

1150  MeV 

a0  =0.  1  34f m 

a,*=0.342fm  a  =0.275fm 

1  s.  0 

a  =0.1 6f m 

ex 

r  =1. 047f m 

0 

r.=1.117fm  r  =0.895fm 

1  s*o 

r  =1.04fm 

ex 

r  =1.36  fm 

r.m.s.  =  1.38  f m 

c 


■ 


' 


- 


49 


Table  3 

Optical  Model  Parameters 


Ep 

(MeV) 

v0 

(MeV) 

Wv 

(MeV) 

Vs.o 

(MeV) 

Ws.  o 

(MeV) 

Vex 

(MeV) 

Jr/A 

M  eV  -  f  m  3 

Ji/A 

MeV-f m3 

100 

6.98 

13.29 

12.71 

o. 

-1.44 

185 

237.1 

156 

6.16 

12.17 

7.88 

-5.09 

.  163 

163 

217.1 

203 

3.45 

10.31 

6.08 

-3.69 

9  1.4 

183.9 

312 

2.22 

14.78 

4.23 

-4.42 

58.87 

263.4 

350 

2.08 

37.81 

3.76 

-3.47 

14.26 

281.3 

560 

-3.65 

50.64 

2.99 

-5.69 

-18.89 

406.4 

650 

-7.34 

64.64 

3.42 

-4.32 

-.  15 

-37.99 

518.7 

720 

-8.38 

80 

2.  52 

-5.24 

-43.38 

642.0 

800 

-11.22 

86 

1.91 

-4.5 

-58 

690.2 

1030 

-19.02 

115 

1.85 

-4.92 

-98.45 

922.8 

1150 

-42 

127.9 

3. 

-5.83 

-217.  4 

1026.5 

CHAPTER  5 


DISCUSSION  AND  CONCLUSION 

The  results  of  the  Optical  Model  analyses  were  detailed 
in  chapter  4.  In  this  chapter  we  will  discuss  the  results 
obtained  and  compare  them  with  the  results  from  other 
methods  of  analysis. 

5.1  Comparison 

Goldstein  et  al.  (  GO70  )  used  the  Optical  Potential  to 
analyze  the  100  MeV  data,  the  same  method  was  employed  by 
Comparat  et  al.  (  C075  )  for  proton  energies  100  and  156 
MeV.  The  latter  also  incorporated  the  147  MeV  polarization 
data  with  the  156  MeV  cross-section  data  to  control  the 
spin-orbit  variations.  The  Optical  Potential  parameters 
obtained  by  these  authors  are  listed  in  table  4.  The  guality 
of  their  fits  obtained  are  comparable  with  our  results  and 
some  general  trends  are  observed. 

1. )  The  exchange  potential  is  necessary  for  the  reproduction 
of  the  large  angle  data. 

2. )  The  sign  of  the  exchange  potential  is  negative  around 
100  MeV. 

3. )  The  r.m.s.  radius  obtained  for  the  exchange  potential  is 
larger  than  the  r.m.s.  radius  for  the  real  central 
potential . 

Comparat  et  al.  IC075  )  obtained  quality  fits  using 

either  a  pure  volume  absorption  potential  or  a  pure  surface 
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absorption  potential  for  the  156  MeV  data.  This  result 
indicates  the  Optical  Model  has  no  preference  between  the 
surface  absorption  or  the  volume  absorption  potential  in 
this  energy  range. 

Dymartz  et  al.  {  DY77  )  analyzed  the  data  at  350,  440, 
650  and  1050  MeV  by  an  ’extended*  Glauber  Model.  The  theory 
takes  advantage. of  the  fact  that  the  Optical  Potential  and 
the  Glauber  profile  function  can  be  related  through  Bessel 
and  Abel  transforms.  The  scattering  amplitude  may  be 
expressed  in  the  form  of  the  impact  parameter  in  the  Fourier 
-  Bessel  integral  which  is  given  as: 


F(q,p)  = 


12. 

4tt 


d2b 


5.1 


where  p,  g  are  the  incident  c. m.  momentum  and  the  momentum 
transfer,  respectively.  In  the  Eikonal  Approximation  the 
scattering  profile  function  nS)  can  be  simply  related  to 
the  interaction  potential  V  (r) : 


r(S)  =  1  -  exp[  —  V(6,Z)  dZ]  . 


For  a  spherically  symmetric  potential. 


r  dZ  r'(/r2+Z2) 

•  1  -  r(/r2+Z2) 


5.2 


V(r)  = 


iv 

7T 


5.3 
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However ,  for  elastic  scattering  of  an  elementary  particle  on 
the  composite  nucleus  in  the  Glauber  model  (  GL59  )  r (b)  is 
given  by 

A 

r(b)  =  i  -  <v|  n  [  i  -  y . (B’-'s . ) ] | v>  5.4 

j=l  J  J 

where  y.  is  the  individual  profile  function  for  target 

J 

nucleon  j,  Sj  being  their  position  in  the  plane  of  impact 
parameters  and  ¥  is  the  nuclear  target  ground  state  wave 
function.  They  evaluate  the  profile  function  described  in 
Eg.  (5.4)  and  work  backward  to  obtain  a  corresponding 
cptical  potential.  The  potential  obtained  is  then  inserted 
into  the  Schrodinger  Equation  to  obtain  the  'exact* 
scattering  matrix.  Using  a  Gaussian  form  factor  for  4He  and 
a  spin  and  iso-spin  independent  profile  function  these 
authors  reproduce  the  forward  angle  data  extremely  well  for 
proton  energies  at  350,  650,  and  1050  MeV,  However,  only  a 
qualitative  fit  was  obtained  for  440  MeV  data.  The  quality 
of  their  results  agrees  very  well  with  cur  results  at  the 
forward  angles.  Since  they  did  not  include  any  exchange 
effects  in  their  analysis  they  failed  to  reproduce  the  large 
angle  data  completely.  These  authors  also  compared  their 
results  with  the  standard  Glauber  Model  calculations  and 
found  that  the  latter  produces  too  deep  a  minimum. 

Kujawski  (  KU7 0  )  used  the  first  and  second  terms  from 
the  Optical  Model  expansion  formulated  by  K.M.T.  (  KE59  )  in 
the  Schrodinger  Equation.  He  combined  the  587  MeV 


. 


' 


■ 


53 


cross-section  and  540  MeV  polarization  data  for  his  analysis 
at  630  MeV.  He  also  analyzed  the  1  GeV  data  from  Palevsky  et 
al.  {  PA67  ) .  Including  both  spin  and  iso-spin  effects, 
resonable  fits  were  obtained  for  the  data. 

Clark  et  al.  (  CL73,  AR76,  SA73  )  used  a  purely 
phenomenological  Optical  Potential  with  only  real  and 
imaginary  parts  : 


U  =  (V+iW)f(r) 


5.5 


with  the  form  factor  f (r)  given  by 


5.6 


where  w,  c  and  z  are  parameters  obtained  from  fitting  the 
electron  scattering  experiment  (  FR67  ) .  The  potential  is 
inserted  in  the  Dirac  Equation  as  the  fourth  part  of  a  four 
vector  potential,  with  the  vector  part  of  the  potential 
taken  to  be  zero.  These  authors  neglected  the  spin-orbit 
effects  in  their  analysis.  They  fitted  the  cross-section 
data  at  587,  600,  650,  720,  1050  and  1150  MeV,  and  obtained 
excellent  results. 

Auger  et  al.  (  AU76  )  used  the  extended  Glauber  Model 
by  Wallace  i  WA71  ) ,  the  scattering  amplitude  is  given  by 


F(q )  =  ikR(q)jelc^  [  1  -{1  -  r(S-s)p(r)  d3r} 
xAe1xc^]  d2b 


5.  7 
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where 

R(q )  =  e<r2>^6A 

J  •  o 

is  the  c.m.  correction  and  the  function  Xc(b)  denotes  the 
Coulomb  phase.  They  used  a  central  Gaussian  form  factor  and 
a  spin-dependent  profile  function  given  by 


r(?-s)p(r)  d3r  =  T (b)  +  a- (£xi<)Ts(b) 


,7> 


5.9 


with 


and 


1 


T(b)  =  j£-  JQ(qb)  A(q)  S(q)qdq 
ck 


5.  10 


1 


T  (b)  = 

S  k 

o 


JQ(qb)  C (q )  S(q)qdq 


5. 1 1 

where  A  tg )  and  C  ig)  are  standard  amplitudes  from 
nucleon-nucleon  interactions  and  S(g)  is  the  nuclear  form 
factor.  Subsituting  the  above  profile  function  into  the 
Scattering  amplitude 


F(q)  =  F-j  (q)  +  a*nF2(q)  .  5. 12 

The  cross-section  and  polarization  are  then  given  by:- 


da/dft  =  |  F-j  ( q )  | 2  +  |  F^  (q )  | 2 

an  d 


5.13 


P  =  2Re(F;F2)  /  ^ 


5.  14 


Auger  et  al.  (  AU77  )  included  all  the  available  data  from 
587  to  650  MeV  for  their  analysis  at  600  MeV  and  also 
analyzed  data  at  0.72,  1.15  and  24  GeV.  They  were  able  to 
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reproduce  the  experimental  cross-section  data  and 
polarization  data  up  to  t  =  -  •  2  iGeV/c). 2  •  Beyond  this  region 
the  calculation  failed  to  reproduce  the  shape  of  the 
polarization  data. 

Wallace  et  al.  (  WA71,  77  )  used  the  extended  Glauber 
model  to  analyze  all  the  1  GeV  data  available.  In  their 
calculations,  they  included  the  effect  of  the 
A  -intermediate  state  (  A  is  the  well  known  3-3  resonance 
in  the  it  -N  system,  where  a  fast  nucleon  becomes  a  A  in  one 
collision,  then  the  A  propagates  through  the  nucleus  and 
returns  to  a  nucleon  state  in  a  latter  collision  thereby 
re-entering  the  elastic  channel.  )  in  a  spin  and  isospin 
dependent  model.  Employing  a  realistic  ♦He  form  factor  that 
fitted  the  e+*He  elastic  scattering  data,  Wallace  et  al. 
obtained  excellent  fits  to  both  cross-section  and 
polarization  data. 

Rule  and  Hahn  (  RU75a, b,c,  76  )  made  an  Effective 
Channel  analysis  of  the  P+4He  elastic  scattering  data.  The 
intermediate  energy  proton- nucleus  scattering  problem  is 
formulated  in  terms  of  coupled  equations.  In  this 
formulation  the  effects  of  the  inelastic  processes  are 
represented  by  an  average  inelastic  channel.  The  theory 
incorporates  approximately  the  effects  of  non-locality, 
energy  dependence,  re-scattering  and  absorption  of  all  the 
inelastic  channels.  The  spin-orbit  term  has  not  been 
included  into  their  analysis.  Still,  they  obtained  fair 
agreement  with  the  data. 
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For  the  lew  energy  region  Ep  <  400  MeV  the  results 
obtained  in  this  thesis  are  comparable  to  the  results 
obtained  by  other  authors.  Higher  energy  region  results  by 
Wallace  et  al.  (  WA77  )  seem  to  be  the  best.  They  are  not 
only  able  to  reproduce  the  cross-section  data  but,  more 
important,  they  also  are  able  to  reproduce  the  polarization 
data  for  the  whole  angular  distribution  of  the  available 
data.  The  phenomenological  method  employed  by  Clark  et  al.  ( 
AR76  )  for  proton  energies  between  500  MeV  to  1100  MeV  is 
also  able  to  provide  excellent  fits.  In  some  cases  their 
results  appear  to  be  better  than  ours.  However,  one  should 
be  cautious  in  making  any  direct  comparisons  due  to  the 
smaller  angular  range  (  t  =  -.7(GeV/c)2  )  in  their  analysis. 

5.2  R.M.S.  Radii  and  Volume  Integrals. 

Greenlees  (  GR68,  70  )  and  his  collaborators  performed 
an  extensive  analysis  on  all  the  optical  potential  model 
parameters  obtained  from  different  methods  of  analysis.  They 
concluded  that  the  r.m.s.  radius  and  the  volume  integral  per 
nucleon  of  the  real  central  potential  are  well  defined 
quantities  in  Optical  Model  analysis.  The  volume  integral 
per  nucleon  J/A  is  given  by: 


J/A  =  -  f  V(r)  d3r 


5.  15 


A  J 
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The  r.m.s.  radius  is: 

<r2>=  (i- r2y(r)d3r } */2  5>17 

J  V(r)d3r 

=  [  i-(3rJ\2/3+  7irV)  p/2  5.18 

The  Woods-Saxon  form  factor  has  been  assumed  (Eg.  3.3)  for 
both  calculations.  The  above  two  quantities  have  been 
calculated  for  all  the  cases  considered  in  the  present 
analysis.  The  results  of  these  calculations  are  listed  in 
tables  2  and  3. 

The  <r2>1  /  2  takes  on  three  different  values,  they  are 
3.173  fm  ,  1.773  fm  and  1.38  fm  for  energies  from  100  to 

300  MeV,  350  MeV  and  from  400  MeV  to  1100  MeV,  respectively. 

In  fig.  17,  we  show  the  <  r2  >*/ 2  as  a  function  of  proton 
energy  including  those  from  previous  investigations.  We 
observe  that  there  is  a  considerable  scatter  in  the  r.m.s. 
radius  data  for  Ep  <  500  MeV.  A  general  downward  trend  for 
the  r.m.s.  radius  as  a  function  of  energy  is  noted  when 
comparing  our  results  with  the  values  2.055  fm  obtained  by 
Satchler  et  al.  (  SA67  ),  1.86  fm  obtained  by  Thompson  et 
al.  (  TH70  )  from  their  analyses  of  Ep  from  31  to  55  MeV  and 
the  value  1.52  fm  obtained  by  Clark  et  al.  *  AR76  )  for  Ep 
between  500  to  1100  MeV.  In  the  higher  energy  region,  the 

proton  elastic  scattering  data  should  be  sensitive  to  the 

nuclear  structure  due  to  the  short  wave  length  of  the 
incident  proton  .  Therefore,  the  r.m.s.  radius  obtained  from 


- 


I  a 


. 

. 


•  ■ 

'  ‘  ••  /■:  .M<  •  1  '  ' 


' 


- 

. 


58 


the  Optical  Model  analysis  at  this  energy  range  should  be  a 
good  indication  of  the  nuclear  size.  The  result  determined 
in  the  present  analysis  i  Ep  >  400  MeV  )  is  1.38  fra  which, 
is  close  to  the  average  value  1.52  ±  .06  fm  obtained  by 
Clark  et  al.  (  AR76  ) .  It  also  agrees  well  with  the  value 
1.42  fin  for  the  proton  matter  distribution  obtained  by 
Frosch  et  al.  (  FR67  )  and  1.677  ±  0.022  fm  obtained  by 
Sick  et  al.  (  SI76  )  •  During  the  course  of  fitting  the  large 
angle  data  we  found  an  interesting  feature  about  the  r-m.s. 
radius  of  the  exchange  potential.  The  r.m.s.  radius  of  the 
exchange  potential  is  larger  than  the  r.m.s.  radius  for  the 
central  potential.  This  feature  has  been  observed  by  some 
other  authors  at  lower  energies.  As  pointed  out  by  Volta  et 
al .  (  VO 74  )  this  result  is  in  contrast  to  the  One-Channel 
Resonating  Group  Theory  which  predicts  a  strong  and  shorter 
ranged  potential.  The  Exchange  Potential  used,  however,  is 
an  over  simplification  and  hence  a  significant  difference  in 
detail  should  be  expected. 

Passatore  (  PA67,  68,  75  )  used  the  dispersion  relation 
developed  by  Feshbach  i  FE58  )  to  analyze  the  results  from 
phenomenological  optical  model  analyses,  to  develop  a 
criterion  for  showing  how  much  the  hypothesis  of  locality 
affects  the  energy  dependence  of  the  phenomenological 
potentials.  He  concluded  that  : 

1.)  The  real  part  of  the  empirical  potential  depth  V  (E)  is 
monotonically  decreasing  with  energy  and  it  changes  sign 


around  300  MeV 
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2.)  The  imaginary  part  of  the  potential  depth  W.  (E)  is 
monoton ica liy  increasing  with  energy  and  goes  asymptotically 
to  a  constant  value  around  1  GeV. 

Van  Oers  et  al.  (  VA74,  73  )  followed  the  above  idea  to 
analyze  proton  elastic  scattering  for  heavier  nuclei  using 
the  Optical  Potential  Model.  He  found  a  logarithmic  energy 
dependence  for  Jr/A  given  by 

Jr/A  -  J0/A  +  B'nEp  .  5.19 

The  same  relation  is  also  observed  by  Clark  et  al.  ( 
CL 73 ,  AR76  )  and  in  the  present  analysis.  In  fig.  18  we  plot 
Jr/A  against  proton  energy.  The  graph  clearly  indicates  a 
downward  trend  as  a  function  of  energy  and  crosses  Jr/A  =  0 
around  Ep  =  450  MeV.  This  result  is  in  fair  agreement  with 
the  results  obtained  by  Van  Oers'  et  al.  (  VA74  )  from  their 
analyses  of  heavier  nuclei  and  by  Clark  et  al.  (  AR76  )  from 

I 

their  analyses  of  P+*He  elastic  scattering  data  from  500  to 
1000  MeV. 

The  volume  integral  of  the  imaginary  potential  Ji/A 
(fig.  19)  displays  a  temporary  decrease  as  function  of 
energy.  This  phenomenon  only  extends  up  to  312  MeV  and  then 
Ji/A  rises  monotonically  as  a  function  of  energy. 
Considering  the  whole  picture,  the  volume  integral  for  the 
imaginary  potential  does  show  a  general  trend  upward  as 
function  of  proton  energy.  Since  our  results  only  extend  up 
to  1.1  GeV  we  can  not  conclude  upon  the  asymptotic  behaviour 
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°f  Ji/A  as  the  energy  increases, 

Agrawal  et  al,  (  AG74  )  analyzed  a  wide  range  of  nuclei 
from  A  =  4  0  -  208  over  the  energy  range  Ep  =  9.8  -  61.4  MeV. 
They  concluded  that  the  total  (surface  and  volume)  imaginary 
volume  integral  per  nucleon  of  the  Optical  Model  is  a  well 
defined  quantity,  in  the  sense  that,  it  is  independent  of: 

a)  the  form  factor  and  the  geometrical  parameters  employed, 

b)  the  mass  number  of  the  trget  and,  c)  the  energy  range 
cover.  The  same  result  was  also  obtained  from  a  subsequent 
analysis  made  by  Hodgson  (  H076  )  over  a  wider  range  of 
nuclei.  However,  it  is  clear  from  our  result  that  the 
analogous  conclusion  can  not  be  made  in  the  higher  energy 
analysis. 

There  was  a  lot  of  speculation  about  the  necessity  to 
include  the  spin-orbit  potential  for  Ep  >  400  MeV  in  Optical 
Model  analyses,  because  of  its  relative  small  magnitude 
compared  to  the  central  potentials.  We  find  the  strength  of 
the  spin-orbit  potentials  agrees  with  other  authors 
assertions  that  they  are  small  compared  to  the  central 
potentials.  But  the  strength  can  be  enhanced  considerably  by 
the  spin-orbit  coupling  term  (  )  due  to  the  large 
number  of  partial  waves  required  in  this  energy  region.  The 
imaginary  spin-orbit  term  also  found  necessary  in  filling 
the  diffraction  minimum.  Its  presence  also  affects  the  value 
of  the  calculated  reaction  cross-section.  Van  Oers  et  al.  ( 
VA74  )  in  their  analysis  of  heavier  nuclei  found  that  a  sign 
opposite  to  the  sign  for  the  lower  energy  spin-orbit 
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potential  was  necessary  to  fit  the  high  energy  cross-section 
data  which  is  in  contrast  to  our  result.  From  the  present 
analysis  of  polarization  data  the  spin  orbit  term  does  not 
change  sign  over  the  energy  range  covered. 

5. 3  Conclusion 

The  Optical  Potential  Model  with  a  spin-orbit  potential 
does  give  an  adequate  representation  of  P+4He  elastic 
scattering  data  over  the  energy  range  between  100  to  1100 
MeV.  The  includsion  of  an  Exchange  potential  does  improve 
the  fit  to  large-angle  data  (fig.  20)  but  due  to  the  lack  of 
complete  angular  distribution  of  cross-section  data  we  can 
not  obtain  any  information  about  the  energy  dependence  of 
the  Exchange  Potential.  The  r.m.s.  radius  of  the  real 
central  potential  obtained  at  the  high  energy  range  studied 
agrees  with  the  values  obtained  by  other  authors  and  the 
results  from  electron  elastic  scattering  experiments.  The 
volume  integral  of  the  real  and  the  imaginary  potentials 
also  agrees  with  the  conclusions  made  by  other  authors  from 
different  methods  of  analysis  and  the  prediction  by 
Passatore  (  PA67b,  63  )  which  is  based  upon  Feshbach's  ( 
FE58  )  dispersion  theory. 

In  this  thesis  we  established  that  an  Optical  Potential 
with  an  exchange  term  can  provide  a  resonably  good 
reprsentat ion  for  the  ♦He  data  in  the  medium  energy  range. 
The  method  of  fixed  geometeries  seems  to  be  more  convenient 
and  economical  than  other  methods  of  analysis.  The  results 
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obtained  here  are  very  useful  for  D.W.B.A.  calculation  in 
this  energy  range.  It  is  also  possible  to  use  this  type  of 
systematic  analysis  to  predict  good  starting  values  for 
phenomenological  optical  potentials  at  other  energies. 
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Fig  18.  Diagram  illustarting  Jr/A  against  Ep 
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Fig  20.  Diagram  illustrating  the  effects  of  the  exchange  potential. 


156  MeV 


Table  4 


Optical 

Model  Parameters  from 

ot  her 

References 

Tp 

(MeV) 

Vr 

(MeV) 

Wv 

(MeV) 

Ws 

(MeV) 

Vs.  o. 

(MeV) 

Ws.  o 

(MeV) 

Vex 

(MeV) 

Ref. 

100 

11.9 

14.8 

8.  4 

.65 

(G070) 

100 

13.07 

8.  77 

7.53  - 

6.53 

-10.24 

(C075) 

156 

5.58 

15.  16 

9.69 

-7.25 

-3.4 

(CO  7  5) 

156 

6.83 

12,2 

7.72 

-5.53 

0. 

(C075) 

587 

31.31 

*123.77 

• 

(AR77) 

6  50 

22.6 

-123.3 

4. 

(AR77) 

720 

26.  17 

-128.2 

(AR77) 

1050 

47.22 

-126.9 

(AR77) 

1154 

51.31 

-126.54 

(AR77) 

Note  Ref.  (AR77)  has  a  sign  opposite  to  the  convention. 
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APPENDIX 


SEEK 

An  Optical  Potential  Search  Program 
The  automatic  search  for  best  fit  parameters  were 
carried  out  by  the  computing  code  *SEFK»,  which  is  a  well 
known  program  written  a  few  years  ago  by  the  U.C.L.A.  Group 
(SE65) .  The  Fortran  program  is  used  for  carrying  out  Optical 
Model  analyses  of  elastic  scattering  data  of  spin-0  or 
spin-1/2  particles  scattered  on  complex  nuclei.  The  program 
permits  automatic  searchs  over  the  space  of  the  nuclear 
optical  model  parameter  space  yielding  optimum  fits  to  the 
cross-section  and  polarization  data. 

The  program  employs  the  Schrodinger  Equation  to 
calculate  the  elastic  scattering  amplitude  from  which 
cbservables  such  as  crossection,  polarization  and  total 
cross-section  are  obtained.  The  Schrodingrer  equation  for 
the  system  is  given  by: 

(  —  V2  +  v  )y  =  Ey  A.  1 

2u 

where  y  is  the  reduced  mass,  E  is  the  energy  in  the  centre 
of  mass,  and  V  is  the  scattering  potential. 

Due  to  the  necessity  to  fit  the  large-angle  data,  the 
Optical  Model  Potential  described  in  Eq.(A.I)  has  been 
modified  to  include  an  Exchange  Potential.  The  modified 
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potential  has  the  form 

v  =  Vc(r)  +  Vr(r)  +  (-l)\x(r) 

+  1(Wv(r)  +  Ws(  r) )  +  2[Vs#(J(r)  +  iW$.o(r)]a4 

A.  2 


where  r  is  the  relative  distance  and  Vc(r)  r  Vr(r)  ,  Vex(r), 
Wv|r),  Ws(r) ,  Vs.o(r)  and  Ws.O(r)  represents  respectively 
the  Coulomb,  Central  Real,  Exchange,  Central  Imaginary, 
Surface  Imaginary,  Spin-orbit  Real  and  Spin-orbit  Imaginary, 
respectively.  These  terms  are  all  assumed  to  be  spherically 
symmetric.  Separating  the  radial  and  angular  parts  Eqs. 
(A. 1 )  and  (A. 2)  give: 


&U+1 )  +  -i 


Vo 

)^(r)  =  0 


A.  3 


where 

vj  =  vc(r)  +  Vr)  +  (-'|)£vex(r)  +  'twv(r)  +  Ws(r)] 


+  2[V$.0(r)  +  iWs.0(r)] 


J  ^  l 

T-£-l' 


A.  4 


where  P=kr.  and  the  factor  1  correspondes  to  j=l+1/2  while 
-1-1  corresponds  to  j=l-1/2.  At  a  large  distance  r=rm  the 
wave  function  is  matched  to  a  linear  combination  of 
regular  and  out  going  coulomb  wave  functions. 


¥ 


£(rm>  =  F^>krn)  +  C*CS(n.krm)  +  i^tn.krj]  . 


A  •  5 


•'  V  -S'  t  1  "f 
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for  a  spin-1/2  incident  particle,  the  two  independent 
scattering  amplitudes  are: 


00 

A(e)  -  fc(e)  +  l  l  e2io*  [(t+ncf1  +£CJ-1/2]P£(COS0) 

1 


A. 6 


B(e)  =  T-f.  e2la£[C*+1/2  -  cj'1/2]pj(cos0) 


£=0 


A.  7 


where 


fc(e) 


-n  exp [i n  ln(sin20/2  +  2 i a 0 ) ] 

2ksin20/2 


A. 8 


The  differential  cross-section  and  the  polarization  for  an 
unpolarized  incident  beam  is  given  by: 

cr(0)  =  | A ( G)  |  2  +  |  B(0)  |  2  A.  9 

P(6)  =  (A*B  +  B*A)/a(0)  .  A. 10 


The  reaction  cross-section  a  is  given  by: 

r 


—  I  {(*+!)[  Imcf1/2 

k2  l=0  1 


|cf1/2|2] 


+  [  ImC 


£-1/2 

£ 


c£-1/2| } 


A.  11 
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Thus  the  necessary  input  parameters  are: 

Geometrical  Parameters, 

Bcr  Fr,  Ri,  Rs.o,  Rex,  Ar,  Ai,  As.o,  Aex. 

Dynamical  parameters, 

V,  Wv,  Ws,  Vs.or  Ws.o,  Vex. 

The  program  searches  for  a  set  of  best  fit  parameters 
by  minimizing  the  total  x2  defined  by: 


x2(0<) 


X*(8*)  +  X2 (0i ) 
0  P  1 


where 


and 


<{01}  '  lm'  -  P>C)'2 


A.  12 


A,  13 


A.  14 


The  program  is  capable  of  searching  on  all  variables 
simultaneously,:  however,  this  would  led  to  large  comsumption 
of  computing  time  and  requires  a  large  memory  space  from  the 
computer,  therefore  the  number  of  search  variables  is 
limited  to  two  at  any  given  search. 
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